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We study the motion of Brownian particle in modulated media in the strong damping limit by 
using toy model, with special emphasis on the transition from localise to diffusive behavior. By using 
model potential we have seen the localised behavior when the number of minima of the potential 
is finite in the asymptotic time limit. In the limit of infinite number of minima we have seen the 
diffusive behavior. We calculate exactly the diffusion coefficient in periodic field of force. We have 
also studied the transport in commensurate and incommensurate media. 

PACS numbers: 05.90.+m,05.10.Gg 



I. INTRODUCTION 

Thermally activated barrier crossing has been subject 
of research for many decades since the pioneering work 
of Kramers [l| on the subject. A fair amount of atten- 
tion has recently been devoted to the study of complex 
non-equiUbrium systems. These include the case of dif- 
fusion over a barrier in the presence of harmonic force 
and the diffusion over a fluctuating barrier The 
hallmark of the former situation is the phenomenon of 
stochastic resonance. The problem of surmounting po- 
tential barriers 0, IBl has gained importance in other field 
of science such as evolutionary computations and global 
optimization as well. In recent work([6j and extensive 
references therein ), we have studied the barrier crossing 
of a time dependent potential which adiabatically evolves 
in time that sheds some light on the global optimization 
problem. 

The transport and diffusion properties of Brownian 
particles continues to attract enormous interest and ac- 
tivity even though a century has passed since the ap- 
pearance of the famous seminal work of Einstein on 
the subject 0- The motion of atoms, vacancies, exci- 
tations, molecules, molecular clusters and colloidal par- 
ticles on surfaces in an active area of research due to 
its theoretical interest "s*] and modern technological rel- 
evance involving self-assembled molecular film growth, 
catalysis and surface bound nanostructuresQ. More re- 
cently, the problem of modeling molecular motors fl^, i.e., 
microscopic objects moving unidirectionally along peri- 
odic structures, has renewed the interest in the field and 
stimulated much theoretical work devoted to the study 
of the directed motion in a fluctuating environment in 
the absence of bias forces. One of the more recent re- 
search foci concerns the experimental and numerical ob- 
servation that even for large clusters of molecules, long 
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jumps spanning many lattice sites may in some cases 
be the dominant contributor to the motion [ll|. In re- 
cent work, the long, even Le vy-like, motions can be de- 
scribed by ordinary Langevin dynamics in the low fric- 
tion regime. Another related work has dealt with the 
transport properties of particles in a symmetric peri- 
odic potential subject to thermal effects and external 
forces involving important physical applications that in- 
clude Josephson junctions, superionic conductors, col- 
loidal spheres and polymers diffusing at interfaces among 
many others [T3|. Several different groups have studied 
the diffusion of Brownian particle in a periodic field of 
force 13J. However, these studies focus on the fact that 
the motion in periodic potentials is necessarily diffusive 
on very long time scales. Restricting ourselves to the one 
dimensional case, the mean-square displacement is given 
by < [x{t) — x{Q)Y >= 2Dt for large times t. Where D is 
the diffusion constant. One can map this problem onto a 
quantum system, representing the motion of a quantum 
particle in a modulated or random potential J^]. On this 
basis, the long time diffusive properties, including the dis- 
continuous dependence of the diffusion constant on the 
'wavelength' of the inhomogeneities in quasiperiodic me- 
dia, is related to the leading low frequency behavior of 
the density of states in the associated quantum system. 
There are the connections between the classical diffusion, 
localization, and intermittencv . 

In the present work, we study the transport of a brown- 
ian particle in the field of force derivable from a potential 
with N number of minima with same barrier heights by 
using Fokker-Planck dynamics. We calculate exactly the 
Kramers' time that the time scale to approach to equi- 
librium in this case and study how it depends on number 
of minima. By taking different model potential with N 
number of minima and of same barrier heights, we have 
seen that Kramers time is independent of the shape of 
the potential. We have seen the localised behavior of the 
particle in the asymptotic time limit in this case. Then 
we take the limit of iV ^ oo to get the periodic potential 
and study the diffusive behavior. We calculate exactly 
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the diffusion coefRcient in this case in the asymptotic 
time hmit. We also study the transport in commensu- 
rate and incommensurate media. 



II. LOCALIZATION IN THE POTENTIAL 
WITH N MINIMA 

Let us consider, a brownian particle is moving in an 
external potential V{x). The overdamped motion of the 
position variable obeys the Langevin equation 



dx 



(1) 



Where, T is the inverse function, rjit) is the random 
force with Gaussian distribution and its correlation func- 
tion obeys < 77(^1)77(^2) >= 2e5{ti — ^2) 

The time dependent probability distribution P{x, t) of 
the random force x{t) obeys the Fokker Planck equation 



dP _ a f^pdV ^ ^d^P 
dt dx dx dx^ 



which, with substitution, 

P(x,t) =exp{-^} <t,(x,t) 

reduces to Schrodinger like equation 



at 



dx'^ 



IT 
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The time independent equation reads. 



ax"^ 



4e 
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1/) = Ai/i 



(2) 



(3) 



(4) 



(5) 



where, ■0(a;, — e 4>{x,t). 



1. Kramers' time 

a. Case I: Let us start with a model potential with 
number of minima is two. Then generalize to N number 
of minima and periodic potential. 

Consider, the potential 



V{x) = K,(l--), 0<a;<a 
a 



^0(1 



-a < a; < 



Vb(--l), a<x 
a 

-l^o(- + l), x<-a 



(6) 



Now, the Schrodinger like eq.(l5]) becomes with the 
above form of potential 



^e'^ + ^[5{x + a)-5{x)+5{x-a)]xl^n = -k^n (7) 



where, 



Consider, 



(8) 



V'„(a;) = Aoe-'^^ + Boe'^■^ -a < a; < 

= A2e-^^'""'\ x>a (9) 

Now, apply the boundary condition at a; = 0, a namely 
ijip is continuous. 

iij^p' is ordinarily continuous i.e. it suffers a discontinuity 
proportional to the strength of the potential. 
We get the following relation between the coefhcients A2 
at X > a and(Ao, Bq) at —a < x < . 



ae '"^ (1 + a)e* 



l + a a \ ( Ao 
—a 1 — a j \ Bq 



(10) 



where, a = i'nix) is symmetric arround a; = 

gives the relation 



-DO — Aq 

1 — a 



(11) 



With the help of PH)) and we get the recursion 
relation of k 



(12) 



In order get the first excited state eigenvalue we use 

+ 5 (13) 

lae 

and from cq. (|12p we get 

(14) 

Now, from ([5]), we get the first excited state eigenvalue. 

(15) 



2ae 



In the asymptotic time limit Ai gives the dominant con- 
tribution in the time dependent probability distribution. 
Hence the Kramers' time that the time scale to approach 
to equilibrium is 



To 



2aV 



(16) 
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Now, consider the number of minima is three and in- 
tervening two maxima are of equal height Vq. As in the 
previous case, Schrodinger hke equation becomes 

-e^ + —[6ix + 2a)-Six + a) + Six) (17) 
—S{x — a) + 6{x ~ 2a)]ipn — —k'^'i'n 

with ^-k'^- 
Consider, 

^„ = Aoe-''=' + Boe*^^, -a < x < 

= Aie-'''' + Bie''=' , 0<x<a 
= Aje-'^^^"") + Bae'^'f^^^'^), a<x <2a 
= Ase-'^^^-^^'l x>2a (18) 

I 



As in the previous case, by applying boundary condi- 
tions on ■0 and ^/i' at a; = 0, a, 2a, we get the following re- 
lation between coefficients (^3, 0) at a; > 2a and {Aq, Bq) 
at —a < X < 0. 



As 




-ka ^ka 
— ka /-I I ^,.\^.ka 



(1 + a)e'' 



1 - Q 



1 + a 



Bo 



(19) 



From (|19p and (fTT|) ,we get the following recursion relation 
of fc 



V 2aefc / V 2aek J 



1+1- 



2aek 2a'^e'^k? 



By substituting k = + 6 into eq. (|20|) . we get 
4ae 

Hence, the Kramers' time in this case is 



(20) 



(21) 



2ro 



(22) 



Now, we generalize the potential with TV minima. The 
Schrodinger like equation takes the form 



N 



dx'^ 



-N 



(23) 



The il^n takes the following form 

ij^nix) = Aoe"'''^ + Boe''^, -a < a; < 

= Aie-'=(^-(*-i)''\ a;>(i-l)a (24) 

As in the previous case, by applying boundary condi- 
tions at a: = 0, a, ... , (i — l)a, we get the following relation 



between coefhcients (Ai,0) at a; > (i — l)a and [Aq^Bq] 
at —a < X < 0. 

( t ) = (^) (^1) (^2) (Tl) ■ ■ • (ri) (Li) ^0 ) (25) 

when A'^ is even. There are (A — 2) number of matrices 
in between (i?) and (Li) matrices, and when A^ is odd, 
the relation becomes 

( t ) = (^) (^1) (^2) (Tl) ■ ■ • m (L2) ( ) (26) 



where,(i?) = ( + 

X / (1 + a)e-*'' ae*" \ 
1^1^ V -ae-*" (l-a)e*'^ ^' 
(T \ = ( ~ a)e~'=" -cie''" \ 

(^1) = ( ^ a" lT« )' a'^'i (^2) = ( x-o, )• 

Now, using (jlip. we get the recursion relation of fc from 
([25ll and ((26|l and then using ^ = ^ + (5, we find (5 is of 
the following form 



2(Af- l)ea 
Hence, the Kramers' time reads 

_ 2(A^-l)aV ^^„,, 



(27) 



{N - 1)to 

Afro, when A^ > I 



(28) 
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b. Case II: Now, we study the approach to equihb- 
rium by using the different model potential. Consider 
the potential with minima two is of the form V{x) — 
Vo(x^ — 1)^. To evaluate the Kramers' time, we have to 
evaluate the first excited state eigenvalue of the Hamil- 



tonian of the form 



It can be 



evaluated by calculating the ground state eigenvalue of 
its super symmetric partner Hamiltonian of the form 

+ + -T^)] [B II3- With the help of trial 
wave function we can evaluate the ground state eigen- 
value of the partner Hamiltonian of the order of e~^°. 
Hence the first excited state eigenvalue of the original 
Hamiltonian is of the order of e"^". Now the Kramers' 
time, the time scale to approach to equilibrium is of the 
order of e^" . 

Now, consider the potential with three minima and 
each maxima is of the same height as double well poten- 
tial. In the similar way we can evaluate the Kramers' 
time of the order of 26^° , where each maxima gives the 
same contribution to the Kramers' time. 

As a generalization, in case of potential with 'N' num- 
ber of minima and each of the A'^ — 1 maxima is of the 
same height as double well potential, the Kramers' time 
is of the order of {N - l)e^° . 

Hence the Kramers' time reads 



(29) 



Where, tq ^ is the Kramers' time of the potential 
with minima two. 



The probability distribution in this case is given by 

= Peqix)+Y,anMx)e-^J^ (32) 

n 

Where, we have used eigenvalue A„(fc) ^ k^. n and 
k are the band and Bloch index respectively. Hence the 
approach to equilibrium distribution Peq is power law 
type in t in case of periodic potential. 

III. DIFFUSIVE BEHAVIOR IN THE PERIODIC 
POTENTIAL 

3. Diffusion coefficient 

The Schrodinger like equation becomes 

dx^ a 
where, A„ 



Y,i-r+'S{x-na) 



Tpn = -k Ipn (33) 



4a^e 

The general solution is 



^ p^k{x-a) ^Q^-k(x-a)^ a<x<2a (34) 



2. Localised to diffusive behavior 

The probability distribution P{x,t) is given by 

P{x,t) = P,g+ciMx)e^^'' + --- (30) 

where Ai is the first excited state eigenvalue. 
The diffusion coefficient is given by 



D 





d 


— lim 




2 t—*oo 


dt 




d 


— lim 




2 t~^oo 


dt 



<x^ > (t) 



— lim 

2 t — >oo 





-ciAie 



-Alt 



dxx'^(f)i{x) + 



(31) 



As long as A^ is finite, Ai is constant. The diffusion 
coefficient goes to zero exponentially in the asymptotic 
time limit. Now, as A^ — > oo, the potential passes over to 
periodic potential. We will see in that case, the transition 
from localised to diffusive behaviour as i ^ oo. Clearly, 
as A^ — + oo, the excited states eigenvalues tend to zero. 
The ground state becomes degenerate and band appears 
in the periodic potential. 



According to Bloch's theorem, the wave function in the 
cell < a; < a is 



(35) 



At a:; = 0, "0 is continuous and '0' suffers the disconti- 
nuity proportional to the strength of the potential. Hence 
we get the following relations between coefficients. 



F = e 



-ka 



2ake 



2ake 



and 



G = e''" 



(36) 



(37) 



Similarly,matching ai x — 2a, we get the following rela- 
tions 



F = e 



and 



G 



-2ik' a 



-2ik' a 



(1 - ^)Ae^'' 
2ake 



2ak 



Vo 
2akt 



Ae 



ka 



(1 



2afcf 



)Be 



-ka 



(38) 



(39) 
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From (1211), (ISTD,® and ^ we find the recursion re- 
lation of k 



cos2fc'a + (1 -p^)cosli2fca (40) 

Where k' is the Bloch index and k is the band index, 
is give 
we find 



p is given by By substituting fc = ^ + 5 in (HDl) 



d = - 



k'aVo 



1 



e 2sinh^Vb/2e 



(41) 



From 



/-TT — k : we get Ai = „ ^ — . , ij ,„ . 
Hence, the diffusion coefficient D is given by 



D = 



1 



2e2 sinh^ Vo/2e 



(42) 



Which is independent of period of the potential and 
only depends on the height of the potential. It can be 
written in the following form 



where, 



D 



1 



(43) 



< e 



>= 



<e^'/«><e-i'/«> 

^ J^^ e^/'^dx where, in the range < a; < a, V{x) 
Vo(l — x/a) and in the range —a < a; < 0, V{x) 
Vo{l + x/a). 



4- Transport in commensurate and incommensurate media 



has a period 27rn. We will see how diffusion coefficient 
depends on n. 

Consider, the potential Vi{x) is of the following form 



V,{x) 

a 

-K,(l-^), 

-3T4,(l-f ), 
6a 



= -{m-2)Voil- 



mVoil ), 

ma 



<x <a 
a < X < 2a 
2a < X < 3a 
3a < a; < 4a 



), (m — 2)a < X < (m — l)a 



{m — 2)a 

(m — l)a < x < ma 



(44) 



and V2 (x) is of the following form 



V2{X) = Vo'il - — ), 
na 

- Vo'il + -), 
na 



< X < na 
—na < X < 



(45) 



We consider the potential V{x) ~ + V2(x). 

Where, V2(a; + 2a) = V2{x) and Vi{x + na) = Vi{x). V{x) 



Now, < >< e ^ > takes the following form 



2na I _ IV 



l_e 2^-0' ^ 2^^,^, _ ^^+2^^^, 



1 - e-2Vo7" 



1 - 



2na 1 + IV 



1 _ e2^"' / -2,.+2 y^, 



1 _ e2Vo7i 



a na 



a na 



1 - e-2Vb7" 
1 - e2Vo' 



1 _ e2Vo7i 



46) 



In the large n limit, the potential V{x) still periodic, 
and we find the finite diffusion coefficient. 



When n goes to zero, the potential varies very rapidly. 
Potential becomes quasiperiodic and localization of states 
occur. Particle becomes localised in the asymptotic time 
limit, hence the diffusion coefficient goes to zero. 



IV. CONCLUSION 

In conclusion, we have studied how Kramers' time de- 
pends on the number of minima of the model potential. 
We have seen that it is independent of the shape of the 
potential. The localized behavior in this model potential 
is studied when the number of minima is finite. In this 
case we have seen the diffusion coefficient goes to zero ex- 
ponentially in the asymptotic time limit. In the limit of 
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infinite number of minima in our model potential we have 
seen the transition from localized to diffusive behavior. 
Approach to equilibrium in the case of periodic potential 
is power law type in time. In the periodic field of force 
we have calculated exactly the diffusion coefficient of tlie 
Brownian particle. In case of rational periodic mixing we 
have seen the diffusive behavior of the Brownian parti- 
cle. But in the case of irrational mixing we have seen the 
localized behavior of the Brownian particle. 
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